Entropy increase during physical processes for black holes in Lanczos-Lovelock gravity 
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We study quasi-stationary physical process for black holes within the context of Lanczos-Lovelock 
gravity. We show that the Wald entropy of the stationary black holes in Lanczos-Lovelock gravity 
monotonically increases for quasi-stationary physical processes in which the horizon is perturbed by 
the accretion of positive energy matter and the black hole ultimately settles down to a stationary 
state. This result reinforces the physical interpretation of Wald entropy for Lanczos-Lovelock models 
and takes a step towards proving the analogue of the black hole area increase-theorem in a wider 
class of gravitational theories. 



Pioneering work by Bekenstein [l[ , Hawking Q , Davies 
Q and Unruh [3] in the seventies showed that there is a 
consistent manner in which one can associate thermo- 
dynamical variables with horizons in general relativity 
(GR). This association gave substance to the formal con- 
nection between the laws of black hole dynamics and 
thermodynamics. 

A natural question is whether this analogy is a pecu- 
liar property of GR or a robust feature of any generally 
covariant theory of gravity. Pursuing this line of thought 
Wald and collaborators [H 0] established the equilibrium 
state version of first law for black holes for any arbitrary 
diffcomorphism invariant theory of gravity. Comparing 
the form of a differential identity with the first law of 
thermodynamics, the entropy of the black hole was ex- 
pressed as an integral over a space-like cross section of 
the horizon of a local geometric quantity and is identified 
with the Noether charge of Killing isometry that gener- 
ates the horizon. 

The standard results in the case of general relativity 
concerning the entropy of horizons rely, in one way or 
another, on the fact that the entropy is proportional to 
the horizon area. This proportionality does not hold for 
the Wald entropy in more general theories and, therefore, 
it is quite intriguing that many of the results connecting 
gravitational dynamics to horizon thermodynamics still 
allows a natural generalization to a more general class of 
models. 

Recent work suggests that this connection may indi- 
cate a far deeper truth regarding the nature of gravity 
viz. that it could be an emergent phenomena like, for ex- 
ample, fluid mechanics Studies show that this cor- 
respondence, between gravitational dynamics and hori- 
zon thermodynamics, transcends general relativity and 
holds true for a much wider class of theories called the 
Lanczos-Lovelock models of gravity Q. These are the 
only natural generalization of Einstein's theory to higher 
dimension if we insist that the equations of motion should 
not be of degree higher than two. The Lanczos-Lovelock 



gravity is also free from perturbative ghosts [lfj and ad- 
mits consistent initial value formulation. As a result, 
Lanczos-Lovelock theories can be thought of as a natu- 
ral extension of general relativity in higher dimensions. 
On the other hand, while Lanczos-Lovelock models show 
remarkable structural similarity with Einstein's theory, 
the form of the horizon entropy in the Lanczos-Lovelock 
models is quite complicated and in general entropy is not 
proportional to any simple geometric quantity. 

Implicit in the investigations which uses the Wald 
entropy in these theories is the assumption that the 
entropy associated with a horizon behaves like ordinary 
thermodynamic entropy. But, the equilibrium state ver- 
sion of first law for black holes, established by Wald and 
collaborators Q requires the existence of a stationary 
black hole with regular bifurcation surface. As a result, 
from the equilibrium state version of first law, it is not 
immediately clear whether the Wald entropy always 
increases under physical processes, except for black holes 
in GR, in which the "area theorem" asserts that area of 
a black hole can not decrease in any process provided 
null energy condition holds for the matter fields fllj ]. 
The area theorem, in turn, follows from Raychaudhuri 
equation and crucially depends on the contracted Ein- 
stein's equation R a t,k a k b = 8tt T a i,k a k b where k a is the 
tangent to the horizon. Since the entropy of black holes 
is no longer proportional to area in Lanczos-Lovelock 
models of gravity, there is no obvious assurance that the 
entropy still obeys an increase theorem. As a result, 
the question of validity of the second law of black hole 
thermodynamics for arbitrary theory of gravity remains 
an unresolved issue. Except for the case of /(i?)-gravity 
there is no proof of the analog of Hawking's area 
theorem beyond GR. In the quasi-stationary case, an 
argument for second law valid for all diffcomorphism 
invariant gravity theories was given in [l2j ]; but it is 
based on the assumption that the stationary comparison 
version of the first law implies the physical process 
version for quasi-stationary processes. 
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For the thermodynamic interpretation to be valid, we 
would expect horizon entropy to increase when a black 
hole in the Lanczos-Lovelock model participates in some 
physical process, like, e.g., accretion of matter. Recently, 
a direct proof of the physical process version of first law is 
proposed for Einstein-Gauss-Bonnet (EGB) gravity [l3| 
which establishes that the net change of black hole en- 
tropy during a physical process is positive as long as mat- 
ter satisfies null energy condition. 

In this paper, we investigate this question for general 
Lanczos-Lovelock models and show that during a phys- 
ical process, the Wald entropy of stationary black holes 
in general Lanczos-Lovelock gravity monotonically in- 
creases provided the matter stress energy tensor obeys 
null energy condition. As a result, not only the net 
change of the entropy is positive, but the entropy is in- 
creasing at every cross section of the horizon. Our result, 
therefore provides a crucial step towards (possibly) prov- 
ing an analogue of the area theorem in Lanczos-Lovelock 
models. 

Let us start with a brief review of the properties of sta- 
tionary, non-extremal, Killing horizons. (We adopt the 
metric signature (— , +, +, +, ...) and our sign conventions 
are same as those of [14]].) In a D-dimensional spacetime, 
the event horizon is a null hyper-surface % parametrized 
by an affine parameter A. The vector field k a ~ (d\) a is 
tangent to the horizon and obeys geodesic equation. All 
A = constant slices are space-like and foliate the horizon. 
Any point p on such slices have coordinates {A, x A } where 
x A , (A = 2, ••• ,D) are the coordinates of a point on 
A = slice connected with p by a horizon generator. We 
can construct a basis with the vector fields, {k a ,l a ,e'\} 
where l a is a second null vector such that l a k a = —1. 
The induced metric on any slice is = g a b + 2fc( a Zb) 
and k a j a b = — l a j a b- The change of the induced met- 
ric from one slice to another can be obtained from the 



metric evolution equation 14j , 



C>klab = 2 [(Tab 
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where a a b is the shear and 6 is the expansion of the hori- 
zon. If the event horizon is also a Killing horizon [l5| . 
i.e. the horizon generators are the orbits of a Killing field 
£ a = (d/dv) a , which is null on the horizon, then the sur- 
face gravity k of the horizon is defined as £ a V a £ b = 
For stationary spacctimes with a Killing horizon, both 
the expansion and shear vanish and using Raychaudhuri 
equation and the evolution equation for shear, we obtain 



R ab k a k b = ei\lH R -bod = k a k c ^^R abcd = 0. (2) 

Note that, in order to derive these relationships, we have 
only used the fact that the horizon is stationary Killing 



horizon with zero expansion and shear without any fur- 
ther symmetry 

We would like to consider the situation when a sta- 
tionary black hole is perturbed by a weak matter stress 
energy tensor and ultimately settle down to a stationary 
state in the asymptotic future. Since the black hole is sta- 
tionary in the asymptotic future, the vector field £ a is an 
exact Killing vector at late times. The accretion process 
is assumed to be slow such that all changes of the dynam- 
ical fields are first order in some suitable bookkeeping 
parameter e and that we can neglect all viscous effects. 
More specifically, we assume that, 8 ~ er ab ~ 0(e). 

In GR, a concrete example of such a physical process 
is a black hole of mass M slowly accreting matter for 
a finite time and ultimately settle down to a stationary 
state. Then a linearized version of the Raychaudhuri 
equation gives, 



~ R a bk k 

d\ 



-8irT ab k a k b , 



(3) 



where, we have used Einstein's equation to get the second 
equality. If the matter stress tensor satisfies null energy 
condition, i. e. T ab k a k b > 0, the rate of change of the 
expansion is negative on any slice prior to the asymp- 
totic future. Since the expansion vanishes in the future, 
the generators must have positive expansion during the 
accretion process. As a result, the area is monotonically 
increasing in the physical process. Note that, the result 
is crucially dependent on the field equation. As a result, 
the monotonicity of the horizon area is only valid in case 
of GR. Our aim is to prove a same statement for the Wald 
entropy during a dynamical change of the black holes in 
Lanczos-Lovelock gravity. 

We shall now turn our attention to the features of 
Lanczos-Lovelock gravity. As discussed before, a natu- 
ral generalization of the Einstein-Hilbcrt Lagrangian is 
provided by the Lanczos-Lovelock Lagrangian, which is 
the sum of dimensionally extended Eulcr densities, 



[D-l)/2] 



a C u 

UL m>~- m i 



(4) 



where the a m are arbitrary constants and is the m-th 
order Lanczos-Lovelock term given by, 



, [D-l)/2] 
rD _ _ \ ca 1 b 1 ...a m b m r?c 1 d 1 

16tt ■* 2 m ldl --- Cmdm lhl 

m—0 



■R 



c m d n 



(5) 

where R c ^ b is the D dimensional curvature tensor and 
the generalized alternating tensor 5 " is totally anti- 
symmetric in both set of indices. The Einstein-Hilbert 
Lagrangian is a special case of Eq. ([5]) when m = 1. The 
field equation of Lanczos-Lovelock theory is, / (16tt) + 
at m E( m ) ab = (l/2)T ab where, 



rpi _ sria.ioy---am.Qm pc 

( m )i 167T 2" i + 1 j c idi--- c mdm a i fc i 
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and to > 2. For convenience, we have written the GR 
part (i.e. for m = 1) separately so that the GR limit can 
be easily verified by setting all a m 's to zero. 

Spherically symmetric black hole solutions in Lanczos- 
Lovelock gravity was derived in [3, EH and the Wald 
ent ropy associated with a stationary Killing horizon is 



f 



S = 



where the entropy density 



P V7 dA, 



(7) 



/ [D-l)/2] 

p = 1 + 167rma m (D - 2) L (m -i) | . 

\ m— 2 

The integration is over (D — 2)-dimensional space-like 



cross-section of the horizon and ^ D 2 ^L 



is the in- 



trinsic (m— l)-th Lanczos-Lovelock scalar of the horizon 
cross-section. We would like to prove that this entropy 
always increases when a black hole is perturbed by a weak 
matter stress energy tensor of 0(e) provided the matter 
obeys null energy condition. 



The change in entropy is [121 ] . 



A *=ijC(;l + ^) dx ^ dA - 



We define a quantity as, 

e 



dp 
dX 



(9) 



(10) 



In case of GR, is equal to the expansion parameter 
of the null generators. But, in case of a general gravity 
theory, is the rate of change of the entropy associated 
with a infinitesimal portion of horizon (see 12j for similar 
construction in f(R) gravity). We would like to prove 
that given null energy condition holds, is positive on 
any slice in a physical process. To proceed further, we 
note that the change of the (D — 2)-dimensional scalar 
(£>-2)£ can b e thought of due to the change in the 
intrinsic metric. Then, we can calculate this change by 
using the standard result of variation of Lanczos-Lovelock 



scalar. The variation of 



- 2 )L 



(m-l) 



simply gives the 



equations of motion of (to— l)-th order Lanczos-Lovelock 
term in (D — 2) dimensions. Therefore, for a general 
Lanczos-Lovelock gravity, we can write 



dp 



[D-l)/2] 

m=2 
[D-l)/2] 

167rma m (D ~ 2) ^_ 1} C klab , (11) 

m=2 



where we have ignored a surface term which does not con- 
tribute because the sections of the horizon are compact 



surfaces without boundaries. ^ D ~ 2 ^TZ a b is the generaliza- 
tion of Ricci tensor for (to — 1)-Lanczos-Lovclock gravity 
and is given by [IH , 

(-D-2)-r>a 1 (TO — 1) r 0l 6i...o m 6 m (D-2) j^adx 

'H (m-l) — 2 ™ °bc 1 d 1 ...c m d m n aibi'" 



(D-2) R c 



(12) 



Then using Eq. (fTJ, we obtain, 

[D-l)/2] 

= 9 + 16tt Y a r 

m=2 



(D - 2) ^ (m -D^ 
(D-2) 



(D-2)rr?ab 



^■lL-l) a ab ) +0 ( D 2) i( m _l 



(m-l) 



(13) 



We would like to study the rate of change of along the 
congruence using Raychaudhuri equation and the evolu- 
tion equation of shear |l4(. We are only interested in 
quantities first order in perturbation over a background 
stationary spacetime. Therefore, when we encounter a 
product of two quantities X and Y, to extract the part 
linear in perturbation, we will always express such a 
product as, 



XYnX^Y^+X^Y^ 



(14) 



where X^ is the value of the quantity X evaluated 
on the stationary background and X^ is the per- 
turbed value of X linear in perturbation. Note that, on 
the stationary background, Raychaudhuri equation de- 
mands R l f b ) k a k b = and since T ( ab ] k a k b = 0, we have 
E^ ab k a k = 0. Also, to simplify the calculation, we use 
diffcomorphism freedom to make the null geodesic gen- 
erators of the event horizon of the perturbed black hole 
coincide with the null geodesic generators of the back- 
ground stationary black hole 24|. 

Using the perturbation scheme mentioned above and 
the evolution equation of 9 and a ab to linear order as 
dO/dX « -R {P) k a k b and da ab /dX « C { a P c ] db k c k d and fur- 
ther using the relationships in Eq.Q for the background, 
the evolution equation of to linear order in perturba- 
tion can be written as 

^ = -8nT ab k a k b + V ab k a k b , (15) 



in which we have defined 

[D-X)/2] 

V ab k a k b = J2 16 



na r , 



[rn)ab 



2m (D-2) F (B)ab R (P) , c , d 



(16) 



Here, we have used expression of the perturbed Weyl ten- 
sor in terms of curvature and Ricci tensors and the re- 
lation ^E {m)ab = ^n (m)ab - (l/2) 7af) ( D - 2 )L (m) . 
We will next prove that the first order part of V ab k a k b 
vanishes identically. To show this let us start with the 
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first term in Eq. (|16[) and write its first order perturbed 
part as: 



E 



(P) 

(m)ab 



k a k b 



167T 2" l + 1 J c i d i--- c m-ld m -icd 



R 



(B)cidi 
ai&i 



? (P)cdj 
h ab 



(17) 



Then, we first expand the background curvature tensors 
in the basis {k a ,N a ,^} on the horizon and use Eq.©. 
We also use the fact that due to the antisymmetry of the 
generalized alternating tensor 5", any component of a 
curvature tensor along the direction of the generator of 

(P) h 

the horizon in the expression of E^ ab k a k will not con- 
tribute. These constrains ensure that the only surviving 
contribution will be from the transverse components and 
we will finally obtain, 



E 



(P) 

(rn)ab 



k a k b 



m 1 ciA 1 B l ...ab 
16n 2 m+1 iC\Di...cd 



(D-2) R (B)C 1 D 1 



R { a P b )Cd Vh 



(18) 



where we have the fact that for stationary spacetimes \2l 



la lb Ip I q ^mn 



(D-2) r>(B)cd 
1«1 



(19) 



which holds for any spacelikc cross section of the station- 
ary horizon. Next, we use the technique in 23| to write 
the alternating tensor in a factorized form as, 



°jC 1 D 1 ...C m - 1 D m ^ 1 cd lt ab 



Wk, 



— 4 c°j C 1 D 1 ...C m - 1 D m _ 1 d a ab K K i 

Using this, we finally get, 



E, 



(P) 

(m)ab 



k a k 



9m < - D -^F {B)ab R (P} k c k d 

1771 (m-l) U acbd K h ■ 



(20) 



(21) 



Eq. (|2Tj) immediately shows that the first order part of 
T> a bk a k b vanishes identically and we finally arrive at, 



~d\ 



-8irT ab k a k b + 0(e 2 ). 



(22) 



Eg. (1221) shows that if the null energy condition holds, the 
rate of change of is always negative during a slow clas- 
sical dynamical process (i.e. ignoring the terms which 
are higher order in the perturbation) which perturbs the 
black hole and leads to a new stationary state. Since, 
the final state is assumed to be stationary, both 9 and a 
and as a consequence, vanishes in the asymptotic fu- 
ture. Hence, we can use the same argument as with the 
expansion parameter in case of GR to conclude that 
must be positive at every slice during the physical pro- 
cess. As a result, we conclude that the horizon entropy 
of black holes in Lanczos-Lovelock gravity is a mono- 
tonically increasing function during any quasi-stationary 
physical process, i.e. 



dS 
dX 



> 0. 



which is what we set out to prove. 

In case of a dynamical scenario, it is possible to write 
down several candidates for the black hole entropy be- 
yond GR [2fjj], such that all the expressions have same 
stationary limit. We have actually chosen a particular ex- 
pression and the validity of Eq. (|2"5|) favors such a choice. 
In fact, in ref. a local and geometrical prescription 
for the entropy of dynamical black holes is proposed. 
This proposal is based on a boost invariant construc- 
tion and agrees with the Wald's Noether charge formula 
for stationary black holes and their perturbations. In- 
terestingly, for Lanczos-Lovelock gravity, the entropy ex- 
pression used in this work matches with the expression 
obtained from the boost invariant construction. Conse- 
quently, our result provides a strong justification in favor 
of the prescription for dynamical entropy as proposed in 
ref. Q. This may also be important to decide the right 
candidate for the entropy of non stationary black holes 
for non Lanczos-Lovelock gravity models. 

From quantum gravitational point of view, a natural 
interpretation of the- black hole entropy is that it counts 
the micro-states of the black hole. For any reasonable 
theory of gravity, which has stable black holes, the den- 
sity of states must be of the form exp(SV), where Sw 
is the corresponding Wald entropy. For example, in the 
context of string theory, it has been shown at least for 
extremal and near-extremal black holes, that the micro- 
scopic computations exactly matches with the Wald for- 
mula 
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28j . Hence, it is quite desirable that the Wald 
entropy satisfies an increase theorem. 

Some obvious further investigations suggested by this 
work are the following: First, one would like to relax 
the quasi-stationarity physical process assumption and 
calculate the full change of the Wald entropy along the 
horizon to understand the validity of classical second law 
for Lanczos-Lovelock models. The possible conclusions 
depend crucially on the signs of the higher order terms in 
Eq. ([22|) . As in case of GR, if all the higher order terms 
are negative, this would imply that has to decrease 
monotonically. Further, it cannot be negative on any 
cross section of the horizon which otherwise will lead to 
—> — oo and hence the existence of a caustic on the 
event horizon which is prohibited by the fact that the 
event horizon is future complete [2!| . Then, to avoid the 
contradiction, must be positive on any arbitrary slice 
of the horizon which would lead to classical second law 
for Lanczos-Lovelock gravity. 

The second issue worth exploring is the following: A 
crucial assumption in our derivation is that there exists 
a quasi-stationary physical process in which the black 
hole ultimately settle down to a final stationary state. 
Although such an assumption is quite reasonable, one 
should not overlook an extreme possibility that the black 
holes in a general Lanczos-Lovelock gravity may not be 
stable under a small perturbation. While the linear 

(23) 

stability around flat spacctime of the Einstein-Gauss- 
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Bonnet gravity is demonstrated in [30|, the positive en- 
ergy theorem has not been extended to a general Love- 
lock theory and even if such an extension is possible, 
there may be other instabilities. This requires further 
investigation. 

Finally, we would like to note that the techniques used 
in this work are specific to Lanczos-Lovelock gravity. As 
a result, it would be worthwhile to find a general ap- 
proach which can answer whether classical second law 
holds in a physical process for any diffeomorphism invari- 
ant gravity theory or applies to a special class of action 
functionals. 
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